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Abstract 

The nature of the interaction of a soliton with an attractive well 
is elucidated using a model of two interacting point particles. The 
system shows the existence of trapped states at positive kinetic en- 
ergy, as well as reflection by an attractive impurity, as found when a 
topological soliton scatters off an attractive well. 
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1 Introduction 



Topological solitons arise as nontrivial solutions in field theories with non- 
linear interactions. These solutions are stable against dispersion. Topology 
enters through the absolute conservation of a topological charge, or winding 
number. U 

It is for this reason they become so important in the description of phe- 
nomena like, optical self-focusing, magnetic flux in Josephson junctions|| or 
even the very existence of stable elementary particles, such as the skyrmion 
[0, |], as a model of hadrons. 

Interactions of solitons with external agents become extremely impor- 
tant. These interactions allow us to test the validity of such models in real 
situations. 

In a previous work M the interaction of a soliton in one space dimension 
with finite size impurities was investigated. 

In the works of Kivshar et al.|| (see also ref. || [7[]), it was found that the 
soliton displays unique phenomena when it interacts with an external impu- 
rity. The existence of trapped solutions for positive energy or, bound states 
in the continuum, is a very distinctive effect for the soliton in interaction 
with an attractive well. 

We can understand the origin of impurity interactions of a soliton by look- 
ing at the impurity as a nontrivial medium in which the soliton propagates. 
An easy way to visualize these interactions consists in introducing a nontriv- 
ial metric for the relevant spacetime. The metric carries the information of 
the medium characteristics. 

A 1 + 1 dimensional scalar field theory supporting topological solitons in 
flat space, immersed in a backgound determined by the metric g^ v in a min- 
imal coupling to the metric, is given by 



g^-d^d^-U^)} (1) 
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where g is the of the determinant of the metric, and U is the self-interaction 
that enables the existence of the soliton. For a weak potential we have||] 



#00 ~ 1 + V(x) 
9u = -1 

9-n = 9i-i = (2) 

Where V(x) is the external space dependent potential. The equation of 
motion of the soliton in the background space becomes 



d 2 (j) r -^Sl_\ r-^_ 



+ 9oo^- = 0. (3) 



dt 2 dx dx_ 

This equation is identical, for slowly varying potentials, to the equation of 
motion of a soliton interacting with an impurity V(x). Impurity interactions 
are therefore acceptable couplings of a soliton to an external potential. It 
is also the only way to couple the soliton without spoiling the topological 
boundary conditions. The source term generated by the metric is essentially 
a space dependent mass term. 

The interaction of a soliton with an attractive impurity shows, however, 
some puzzling effects^, A soliton can be trapped in it, when it impinges 
onto the well with positive kinetic energy. Energy conservation demands 
that the soliton fluctuates and distorts in trapped states inside the well. 
Even more counterintuitive is the fact that the soliton can be reflected by 
the well. 

Neither of these effects are possible for classical point particles. The 
difference must obviously be due to the extended character of the soliton. 

This was indeed put in evidence in the works of Kivshar et al.||. It was 
shown there that the bulk soliton behavior may be reproduced qualitatively, 
by taking the center of the soliton as a time dependent collective coordinate 
coupled to the major excitation modes of the soliton in the well, the impurity 
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mode, and a shape distortion mode. The first mode is excited inside the well 
only and appears as an oscillating packet centered at the well. The shape 
mode is an excitation that accompanies the soliton with a time dependent 
amplitude along its scattering, and accounts for the distortions of the free 
soliton. The dynamics of the soliton was then replaced by the equations of 
motion of three classical particle-like excitations: the center of the soliton, the 
amplitude of the impurity mode and the amplitude of the shape mode. For 
a S function type of well, it was found that the system behaves analogously 
to the soliton. Kinetic energy of the soliton center can be transferred to 
the impurity and the shape modes. The system of three effective degrees 
of freedom can resonate inside the well and be trapped. Reflection by the 
attractive well is also observed. The system captures the essential features 
of the behavior of the soliton. 

However, both the interactions between the collective degrees of freedom 
and the dynamics of the system are derived from the soliton itself. Moreover, 
the behavior was shown to hold for 5 type of wells only. 

Later it was shown M that all the features of the scattering of a soliton 
show up in the case of a fine size well too. 

In the present work we show that the effects are not specific to the soliton 
generated dynamics and interactions. 

We are able to reproduce here all the abovementioned effects with a clas- 
sical model for an extended object, regardless of the soliton dynamics. In 
the next section we show that a simple classical model exhibits the same be- 
havior as the soliton. It can be trapped and reflected by an attractive well. 
The system will also show chaotic behavior. 

The classical system can serve as a nice introductory example for the 
surprising behavior of solitons. It demands only basic lagrangian mechanics 
knowledge, but, it has many interesting features that are easily visualized by 
undergraduate students. It may serve also as an example of chaotic behavior 
in classical mechanics. 
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2 A classical model of trapping 



Kinks and other topological deflects are sometimes used to model classical 
systems of masses. We will here find, that the analogy is more than math- 
ematical. The very behavior of the soliton is exactly the one observed on a 
two-body system interacting with an external potential. 

In ref . || , the dynamics of a soliton interacting with a well was also studied 
by selecting the collective coordinates of the soliton center, its shape-mode 
excitation amplitude and an impurity mode. Using this scheme, it was found 
that this classical system shows an analogous behavior to the soliton itself. 
Namely, trapped states and reflection by the attractive impurity. The spe- 
cific dynamics of the soliton in terms of its collective coordinate, and its 
excitations including the shape mode were crucial for the obtention of the 
abovementioned results. 

We show here that the effects are generic, any system resembling the 
behavior of the soliton, mainly, its extended character, will indeed yield sim- 
ilar results. This is an unexpected situation whose motivation arose entirely 
from the behavior of the soliton and its relevance goes beyond it, as will be 
explained below. 

In order to justify the connection with the behavior of the model we 
appeal to a physical scenario in which both Sine-Gordon solitons and Kinks 
pL0|| Sine-Gordon solitons arise in large Josephson junctions and in the 



arise. 



motion of dislocations in a one-dimensional crystal. Kinks arise in the latter 
also when the substrate potential is nonperiodic. 

Consider the Hamiltonian of dislocations in a crystal with nearest neigh- 



bor interactions 10 



1 ,du n 2 G . , 2 , . 

2 m ^~dt' ~2 ^ Un+1 ~ Un ' ^ Un ' 



(4) 



Where u n is the displacement of the dislocation at site n, G is the spring 
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constant between particles and V is a site potential generated by the sub- 
strate chain of fixed particles upon which the mobile dislocations move. 

The above Hamiltonian supports solitons in the continuum, strong cou- 
pling limit. In particular for kinks, only a few dislocation centers are needed 
to generate the desired effect of the moving soliton. The minimal set would 
then be a couple of dislocation 'particles' moving along the substrate. Now 
suppose the above model is applied to a substrate for which the parameters 
of the substrate potential vary. This is analogous to the variation of the 
metric in the description of the previous section. In such a scenario we have 
to modify the substrate potential by adding a local interaction at fixed sites 
in the lattice. This is essentially the procedure in soliton-impurity scatter- 
ing. If the effects found by Kivshar et al.@, are indeed based on the above 
simple dislocation model, then these should appear clearly when a couple of 
dislocation degrees of freedom scatter off an external potential. We will see 
below that this borne out in a simple two-particle model that imitates the 
soliton behavior. 

Consider a system of two classical point particles connected by a massless 
spring, two of the dislocations above, and a repulsive force between them 
needed to prevent their collapse to zero size that subsums the behavior of 
the rest of the chain of dislocations. With only two degrees of freedom, we 
are eliminating the rigidity of the chain, thereby introducing the spurious 
possibility of complete overlap between the two sites, which does not occur 
when the chain is infinite. Hence the need for a repulsive interaction. 

The above simplistic model is not directly related to the collective coordi- 
nate treatment of ref. || deliberately. The aim is to show that the particular 
dynamics of the soliton is not the cause of the peculiar phenomena previously 
found. 

When each particle in the system is allowed to interact with an external 
potential, we are imitating the impurity force or the local change in the 
dislocation potential. 
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The classical nonrelativistic one- dimensional lagrangian for the system of 
equal masses mi = wi2 = 1 becomes: 

tsys = f + f " * _ + V{xi) + V(X2) (5 ) 

For the potential well we take 

V(x) =Ae~P x2 (6) 

Although any finite size well may serve for this purpose. We prepare the 
two-particle system at rest at a large distance far away from the well with 
an initial speed v. The equilibrium interparticle separation is ro n+2 = 
We here use n = 2. 

The equations of motion are not solvable analytically. However we can 
show that for a well large compared to the equilibrium distance the system 
may be trapped and oscillate inside it. Transforming to relative and center of 
mass coordinates, r = X2 ~ X1 , R = Xl + X2 and using the ansatz r = r + 5(t), 
with 5 a small parameter, we find the equations of motion near the center of 
the well R = 



8 +2 k 5 + 2 A e~ r ° 13 (3 (r + 5) = 

R + 2APR e- r ° = (7) 

Where we have used (3 r§ << 1, a wide well as compared to the equilib- 
rium distance of the system. Passing to a new coordinate 



1 + -h e 

the first of equations [7| becomes 

e + 2 k e + 2 A e" r o P f3e = (9) 
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It is clear from the above equations that the center of mass coordinate 
of the system oscillates around the center of the well, while the relative 
coordinate oscillates too with the small amplitude e. Moreover the system 
shrinks inside the well. The oscillations of the center of mass coordinate r 
compensate for the loss of kinetic energy of the system that impinged from 
infinity with a fixed relative separation. The above treatment demonstrates, 
that at least there is room for the trapping to occur. 

We now proceed to show the numerical results of the exact calculation of 
the development of the model. 

Using the numerical method used in ref. || we can find the outcome of 
the scattering events as a function of the initial speed. The system starts 
with a certain initial center of mass velocity v far away from the well. The 
system is prepared with the relative separation of equilibrium tq and the 
outcome of the scattering is monitored as a function of the initial conditions. 

Figure 1 exemplifies the results for the choice of parameters k = 1, a = 
1, n = 2, A = 2, (3=1. 

Quite unexpectedly, it is found that the system behaves exactly like the 
soliton. 

The system can be trapped Vf ina i = 0, reflected, Vf ina i < or transmit- 
ted, v fi na i > through the well by varying the initial speed. 

When the system is trapped, it oscillates with a null average speed, the 
kinetic energy stored in the vibrational and deformation modes. 

Minute changes of the initial speed around a value leading to a trapped 
state, may generate reflection or transmission events. 

The effects are independent of the functional dependence of the interac- 
tions and external potential, as well as the values of the parameters. 

In figure 1 we used a grid for v of dv = .001. Using a finer grid, each region 
of reflection-transmission unfolds to more islands of trapping, reflection and 
transmission. 

Finer and finer grids show more and more structure. 
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Figure 1: Final velocity of the two-particle system as a function of the initial 
velocity for the parameters k — 1, a — 1, n — 2, A = 2, f3 = \ with a 
velocity grid dv=.001 



Figure 2 shows a detailed expansion of the velocity range around v=.12 
with a grid dv=.0002. The system is chaotic, an infinitesimal change in the 
initial speed produces diverging results. 

Many of the phenomena related to chaotic behavior may be identified in 
the system, namely scaling, bifurcation and perhaps even fractal structure. 

It is now safer to claim that the unexpected behavior of a soliton inter- 
acting with an attractive well may be traced back to its extended nature. If 
we regard each <p(x) as a classical pointlike object we will find interactions 
between neighboring particles of attractive and repulsive character. The ba- 
sic attractive interaction is provided by the space derivative of the soliton 
lagrangian and a piece of the self-interaction potential, while the repulsive 
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Figure 2: Same as figure 1, but with finer velocity grid dv = .0002 

interaction is provided by the latter and the coupling to the remainder of the 
soliton or linear chain ina discrete model. 

3 Final remarks 

The simplest implementation of the system studied would be a toy-like sys- 
tem of two masses tied-up to a spring sliding on a frictionless table with a 
carved well on it. Two atoms in a molecule scattering off an external Van- 
der-Waals potential might show the same effects in a quasi-classical approxi- 
mation. However, quantum effects can blurr the picture due to interference. 

Turning the process backwards: an extended object, may it be a soliton 
or a classical assembly of bound particles, in a trapped state, can suddenly 
be freed from it provided some random interaction causes the reversal of the 
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process of trapping, a process reminiscent of the decay of metastable states 
in quantum mechanics. The concept of trapping in general is not discussed 
in the classical mechanics literature. Soliton physics has taught us that an 
extended object made of particles in interaction has a much richer behavior 
than expected. 
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